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SUMMARY 

The damping— in-roll parameter Clp is calculated theoretically 
for triangular wings on cylindrical bodies and for a class of wings 
with swept-back plan forms. The analysis is based on the usual assump- 
tions of linearized compressible— flow theory together with the added 
restrictions that at the free— stream Mach number Mq the product of 
1-Mo 2 and the streamwise velocity gradient is small.- The accuracy of 
the results is tested by a comparison with the exact solution of the 
linearized equation for a triangular wing and for such plan forms is 

shown to increase as reduced aspect ratio (y |1-M 0 2 | times aspect 
ratio) decreases. 


INTRODUCTION 

Theoretical advancements in the study of load distributions over 
three-dimensional wings in compressible— flow fields have been achieved 
almost entirely under assumptions leading to the linearization of the 
basic differential equations of flow. In this manner the analysis is 
resolved in general into the problem of determining solutions of elliptic . 
and hyperbolic partial differential equations in three dimensions, the. 
character of the equations being fixed by the magnitude of the free— stream 
Mach number "Mo. It is possible, however, to study also a class of prob- 
lems associated with the parabolic form of the potential equation; namely, 
cases for which the product of 1-M 0 2 and the chordwise perturbation- 
velocity gradient is small in comparison with the gradients in their 
respective directions of the perturbation velocities perpendicular to the 
chord. • Such a method was developed in reference 1 for small values of 
aspect ratio and was used in reference 2 for swept-back lifting surfaces 
at sonic speeds. In the present report the method will be referred to 
as a slender— wing theory wherein the term slender infers that the ratio 
of the reduced 1 span to the over-all length is small. Hence the results 

^-Reduced span is defined as P times the span, where p is equal to 

J | 1-Mq 2 | . In general, the term reduced implies that any parameter 
which it modifies is multiplied by the factor P. • 
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presented will apply most closely when the wing is long in terms of its 
span, or when the free— stream Mach number is close to unity and the wing 
and body do not violate the assumptions made . 

Previous investigations (references 3, 5, and. 6) corresponding 

to this method of approach include, respectively, a complete analysis of 
all the stability derivatives of a low-aspect— ratio triangular wing, a 
lifting triangular wing with an arbitrary body of revolution, a lifting 
triangular cruciform combination on an arbitrary body of revolution, and 
a lifting swept-back constant-chord wing. In order to assess the accu- 
racy of the results presented in these reports, it is necessary to com- 
pare their approximate solutions with the exact solutions of the linear- 
ized equation. For the case of the rolling triangular wing, this can 
easily be done since the exact solution has been derived in reference 7. 
The results of such a comparison are presented in figure 1 which indi- 
cates that in this case the value of PCi p for the approximate theory 
can be useful up to reduced aspect ratios as high as three. 

The object of the present report is three-fold: First, to find the 

effect on the damping-in-roll parameter Cj p of adding a body of revolu- 
tion to a triangular wing in order to extend the knowledge of wing-body 
interference into the field of lateral stability; second, to find Cj p 
for a particular swept-back wing plan form; and third, to show how the 
damping due to roll can be found for swept wings with arbitrary trailing 
edges or how trailing edges can be calculated from prescribed span load 
distributions. The examples to be given provide sufficient details to 
indicate how other cases can be calculated numerically. 

A list of important symbols is given in Appendix A. 

RESUME OF THE METHOD 

Under the assumption that (l-Mo g )j>xx is small as compared to $yy 
and $> zz , the equation for compressible fluid motion, either subsonic 
or supersonic, becomes ’ 


®YY + ®ZZ = 0 (la) 

where <5 is the perturbation-velocity potential and X is the free— 
stream flow direction. The resulting differential equation (la) shows 
that although 0 can be a function of X its value is a consequence of 
boundary conditions given along lateral strips. We will seek solutions 
to equation (la) for which first, the perturbation velocities vanish at 
infinity; second, there is no discontinuity in potential except across 
a lifting surface or its trailing vortex sheet; and third, there is no 
discontinuity in the vertical velocity anywhere in the field (i.e., no 
airfoil with thickness). It is possible, by applying the two-dimensional 
form of Green's theorem, to write the solution to equation (la) which 
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satisfies all of the aforementioned conditions grid relates w^ 

(0 Z in the Z = 0 plane) to Avt 0 (the jump in sidewash, $ T in the Z = 0 
plane). Under tie added restriction that no point exists in or on the 

flow field such that " g — - / 0, where e is the radial distance from 

— ^ oe 

the point, this solution can "be written 2 




->B(X) Av t dT x 


-B(X) 


The integral equation (2a) will be solved for the case of rolling 
wings. First, however, it is convenient to introduce nondimens ional 
lengths so that the results can he generalized as far as possible. A 
satisfactory nondimens ional form is obtained by using the equations 


C tan 9 


x = £ 


u = C Ut 

v = C v-t tan 6 

w = C w t tan 9 


T r = C v„ tan 9 


C tan 9 


H = p C 2 tan 2 9 


where 9 is the semi vert ex angle of the wing, p is the rate of roll in 
radians per second, C is the root chord, uj., v t , w^.,and v r are true of 

velocity components, and u, v, w, and v r are the transformed components. 
In addition cp(x,y,z) replaces $(X,Y,Z). 


TRIANGULAR WING WITH BODY 

Consider a wing with an unswept trailing edge mounted on a body of 
revolution as shown in figure 2. Applying the transformations given by 
equations ( 3 ) to equation (la), we find that in every plane perpendicular 
to the x axis it is necessary to satisfy the equation 


2 This equation is developed in more detail in reference 2. 
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subject to the boundary conditions 


o 

1) 

for 

r 

= 

0 < \|r < 2it 

wo = ny 

for 

z 

= o, 

a < |y| < b 

O 

II 

> 

II 

> 

for 

r 

= OB, 

0 < \|r < 2* 


(M 


where a(x) and b(x) are replaced by a and b for convenience in 
notation and where r,t are polar coordinates in the yz plane and p 
is related to the rate of roll by equation (3). These boundary condi- 
tions are not of a type such that they can be substituted directly into 
equation (2a), since one condition specifies a radial component which 
is, moreover, out of the z = 0 plane. However, equation (l) is simply 
Laplace f s equation so that by introducing complex variables it is not 
difficult to transform the wing and body section (fig. 3(a)) to a strip 
along the y axis (fig. 3(b)). Such a mapping function is provided by 
the Joukowski transformation, 3 which, if the subscript i denotes 
conditions in the transformed plane, can -be written: 


li-S+f” (5) 

where § and £ x are given by the expressions 

S =7+iZ l (6) 

li = yi + iz i 

By means of equation (5) the wing in the £ plane (the portion of 
the axis for which a < | y | <b, fig. 3) transforms into the section of 
the yi axis in the . £ x plane for which 2a = a x < | y x | < b x = b + a_ . 
Further, the body in the | plane (the curve satisfying the equation r=a) 
transforms into the part of the y x axis for which 0 ^| y x | <a x = 2a. 

In this manne r the boundary on which the data are specified has been 
transformed so that it lies entirely along the y x axis. It is next 
necessary to inspect what the boundary conditions are in the | x plane. 

From the basic theory underlying the use of complex variables, induced 
velocities in the physical and transformed planes are related by the 
expression 

v — iv - (v x — iw x )-^- (7) 

a-! 


3 See, for example, "The Elements of Aerofoil and Airscrew Theory," 
Cambridge University Press, 19^-3 > hy H. Glauert. 
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from which, since in polar coordinates 


a f ’ l 1 '(?) 008 2 ♦] + 1 (?) 


sin 2 t 


(8) 


it follows i m mediately that 


v = Vi 


1 — f £ ) cos 2 \|r 

JTj 


+ Vi 




sin 2 ijr 


▼r 


= *i [l ~ (j) cos 2 i|rj — Vj. sin 2 i|r 

= £vi cos i|r + w x sin tjr 


>( 9 ) 


Through the use of equation (9), the "boundary conditions given by 
equation (4) for the £ plane can be transformed to the plane with 

the result (since i|r = 0 or it) 


w x = 0 for 0 < yi < ax 


— 7^4 - 


l yl 


= w x = 0 for ri = ®, 0 < f < 2jt 


> < 10 ) 


The problem has now bee cone one of solving Laplace's equation for the bound- 
ary conditions given by equation (10). But this problem can be immediately 
solved, provided we can invert the integral equation (2a), which in the new 
notation, 4 becomes 


Vi(yi) 



Av 1 (y 2 )dy g 

7i-72 


(2b) 


4 

Note that the subscript o (indicating conditions in the Z* = 0 plane) 
has been dropped when in the g plane. This avoids a cumbersome 
notation and should cause no confusion, since in the g x plane only 
conditions on the y x axis are to be considered. 
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Under the condition 



vi(yi)dyi 

*/bi 2 -7i 2 


( 11 ) 


reference 8 gives an inversion to equation (2b) which can be written 


= | J t>i 2 -yi 2 r 1 - 

-bi 


vx(y 2 )dy 2 c Q 

—■ " ■■ ■ " — — 4- ■ '■ 

(yx-yaV i*i 2 — y 2 2 */ b x 2 — y x 2 


( 12 ) 


where c Q is an arbitrary constant. 


An inspection of w x as given by the boundary conditions in equation 
(10) shows that it is an odd function of y x and hence will satisfy equa- 
tion (11). Placing w x in equation (12) and integrating gives 


Avi(yi) 


2 v a x 2 — y x 2 V 2 y 


G 2 


— + arc cos 



vyi 


In 


y X */b x 2 -a x 2 + a x 




yi 


A/bi 


•yi 

2 -a x 2 - a x */b x 2 — y x 2 




■yi 


(13) 


where o = -1 for 0 <|yi(<a x , and a = +1 for a x <|y x |< b x . 


(Ik) 


Now to find A<p x on the z x = 0 plane, the relation 

^yi 

/- ^ y2 
-bx 

in which Av x (y 2 ) as given by equation (13) is used with the result that, 
for a x < y x < b x , 


r>y 1 

Acpx = J Av x (y 2 )dy s 


Acp, = 


~*( J 


1 2 ai 

1 + ^ arc cos ^ 


^yi*/ 


V 2 2 

bx -yx 


iL (y x 2 — a x 2 ) 71 ^i 2 - 5 ! 2 

2« bx Vyx^x 2 


( 15 ) 
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where the constant c Q has been evaluated by the conditions that A<p x 
must vanish when yi = b x . This can be written in terms of the original 
coordinates (fig. 2(b)) and boundary conditions (equation (4)), by 
retransforming in accordance with equation (5). Thus, for a<y<b. 


A<Po 



2ab 


b 2 + a 2 



cosh 


(y^+a 2 ) (b 2 -a g ) 
(y 2 -a 2 ) (b 2 +a 2 ) 



(16) 


Equation (l6) represents the solution, in terms of the velocity poten- 
tial, for a rolling wing with an unswept trailing edge on a body of revolu- 
tion; both a and b are functions of x. The loading coefficient ^2 
can be determined by taking the partial derivative of equation (16) 
with respect to x in order to find Au, and then using the relation' 


Ap _ 2Au fc _ 2 Au 

<1 cv D 


(17) 


However, it is not necessary to find the loading over the wing if interest 
is limited to the determination of the total rolling moment L. If a is 
independent of x (i.e., the body is a circular cylinder), then r = (Acp'w, 
and the equation . ■ L * 


dZ = p V 0 T dY 


where dZ is the increment in span loading, can be used to express 

dL = TdZ = P V 0 Y T dY 


Thus 


or 



L = 2p Y 0 C 2 tan 2 6 





T.E. 


dy 


(18) 


where (Aq> ? ) T>E is the value of A9 0 at the trailing edge or where 
1 in figure 2(b). Substituting equation (16) in equation (ifl\ nro 
can .derive, after some rather involve! integratiL^the" resS t ( ^ °“ 
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L 

pv 0 


-1 . 4 * -n r 1 r (a 2 +s 

= — C tan 0 2 Ji 1 ? 

, « \2 [ s 2 


2 . _2 


arc tan — 
a 


-2j- (s 2 — a 2 ) 6a 2 s 2 +s 4 ^ arc tan |j- 


j^a 4 , 1 a 2 ( 0 2 _ 2\ 2 1 

— +2 ^ i 


(19) 


The damping-inr-roll coefficient C ^ based on the wing area including 
the part of the wing masked by the body can now be written, if both 
sides of the equation are multiplied by P -J\ 1-M 0 2 | , in the form 


Ap 

p % ■ - fi 


(l+R ) arc 


tm k 


+ 2B(1-B 2 )(R 4 -6E 2 +1) arc tan i - 


Tt^ 4 + B 2 (l-R 2 ) 2 j- (20) 

where E = - = - (i.e., the ratio of the body diameter to the wing span) 
s S 

anfl Ar is reduced aspect ratio which is also based on the total wing 
area including that part masked by the body. 


When E equals zero, equation (20 ) represents the damping in roll 
for a wing without body. This value, which will be designated (C Jp)y, 
is given by 


»<c hK’-w (21) 

and agrees with the value given in reference 3* The ratio of equation 
(20) to equation (21) represents, finally, the effect of adding a circu- 
lar cylinder about the' line of symmetry of any pointed, low-reduced-aspect- 
ratio wing with an unswept trailing edge. This effect is shown in figure 4. 
The figure shows that the damping in roll is increased a maximum of about 
4 percent at a value of E equal to 0.28. 

SWEET-BACK WING 

The study of the swept-^back wing will be divided into two parts. 

The first of these sections will be concerned with the development of a 
special plan form which is relatively easy to treat analytically, while the 
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second section will be devoted to the examination of plan forms with 
arbitrary trailing edges or vith arbitrary span'' load distributions. 


Special Plan Form 

Consider a wing with a plan form similar to that shown in figure 5. 
The analysis of such a wing is best presented by considering separately 
the three regions shown in figure 5(b) which correspond, respectively, to 
values of x for which x<l, l<x<s, and s<x. 

Region I (x <;!) .— Since the loading in each region is independent of 
the flow at all points downstream of it, the solution for the pressures 
in region I is the same as that, for a triangular wing. Hence region I 
is simply a special case of the preceding section for which the body 
diameter E equals zero. 

• 

Region II (lS-x^s) .— Although each region is certainly independent 
of subsequent ones, the vorticity of upstream regions must be tr an sported 
downstream through all subsequent transverse planes. This increases 
somewhat the complexity of the problem. 

In order to find the solution for the loading in region II, it is 
again necessary to satisfy equation (l), but this time subject to the 
boundary conditions: 

Wo = My for a< |y |< b 

u = w = 0 for r = co, 0 < \|r <2rt V (22) 

Auq = 0 .for. 0 < |y |<a 

. J 

where 5 a(x) and b(x) are, as in the preceding section, replaced by a 
and b. 

It is apparent that no transformation of the plane is necessary for 
the use of the integral equation (2), since the boundary conditions are 
already specified in the z = 0 plane. Making use of the fact that Acp 
must be an odd function of y, and hence its derivative Av 0 even, 
equation (2b) can be written 



5 Whereas in the wing-body section a(x) represented the radius of the 
body, here a(x) represents the distance from the x axis to the 
trailing edge. Since the two sections are entirely separate, this 
should cause no confusion. 
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where 


and 


n = y 


G(x,tj) = 


Av 0 (x,y) 


( 24 ) 


Before the solution to equation (23) can be obtained, however, the behav- 
ior of the flow in the vortex wake, the region for which -a< y< a, must 
be discussed. 


Since an increase in Acp 0 corresponds to a loading Aq> 0 at any 
point behind the wing must be the same as A f Po at the trailing edge for 
the same y. If the equation of the trailing edge is written in either 
of the forms • 


y = a(x) 
x = a*(y) 


( 25 )' 


then Atp 0 (x,y) equals Acp 0 (a*,y) for x>a*. Taking the partial derivative 
of Acpo(a*,y) with respect to y gives 


dAcp 0 (a*,y) _ dAq>o(a*,y) da* + 
ciy da* ^7 


t 


^Rp(x,y) 

dy 


x= a* 


which can also be written in the form 


^^> 1 } A P< a *>y) + Av 0 (a*,y) (26) 

dy 2 dy q 


Applying the Kutta condition, the value of the load coefficient at the 
trailing edge 4&(a*,y) must be zero. Equation (23) can now be written 

= _ i_ r a 2 G*(Tii)dm _ r b 2 G(x,Tji)di 11 (27) 

^ 2ic J n-ni 2n J-, 

O 

where G*(bi) = Gfa*,!^). 


The value of G*(tu) depends, of course, on the equation of the 
trailing edge. It is convenient at this stage of the analysis to choose 
the value of Acp 0 in the wake (since applying the Kutta condition to 


equation (26) yields G*(tji) 


± f BAcp 0 (x,y) "1 
7 L dy J x=a* ’ 


giving Acp 0 is equivalent 
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to specifying G*(qi)) and find the resulting trailing edge rather than 
to proceed in the more straightforward way of first completely specifying 
the wing plan form. However, if Acp Q is expressed in the form 


— - = Biy + Bay 2 + . . . + Bmy m (28a) 

-H 

it will he shown that both the trailing edge and the final value of Cj p 
can he given in terms of the B n 's, and, therefore, quite general types 
of plan forms can he analyzed. Our purpose is to present in some detail 
the simplest example in this section and to develop the general discussion 
later. 

Figure 6 shows, the variation of Acp 0 in the wake behind a triangu- 
lar wing (equation (15) with a = 0, h = l). A reasonable first choice of 
Acp 0 ' in the wake of a swept— hack wing would seem to he the straight-line 
continuation of the slope at the origin of figure 6. This leads to the 
equations 


B — 1 , B2 — Bg — ... — 0 
so that equation (28a) becomes 

G*(ni) = - h/ */iu 
6 

If G(x,tu) is defined by the relation 

G(x,t|i) = Gi(x,t)i) + H(x) 



(28b) 


(29) 


where H(x) is a function to be determined independent of t^, then 
equation (27) can be written 


am 

(n^ni) N /~rir 


Gi( x > T ii) < 3- T li 


It is now possible to invert equation (30) by means of equation (12), 
wherein a 2 now replaces — bi and b 2 replaces bj.. First we notice 
that the condition imposed by equation (ll) is satisfied if 


6 

The choice of such a definition is, of course, not obvious. Its advis- 
ability depends entirely on the simplification which it introduces 
into the subsequent analysis. 
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_ J±_ 


2 * 





H(x) 

2 


At 1 

•/(b%)(,-a 2 ) 


= 0 


or, upon integrating, that 


where K is the complete elliptic integral of the first kind of modulus 
k = The solution for Gi(x,t)) in the interval a <|y|<b becomes 


Gi(x,tj) = ^ V(b 2 -r| ) (■q—a 2 ) 

' a 2 (tj-tii) •/ (b -^i)(tix— a ) 

The single integral evaluates to give zero, and the order of the double 
integral can he reversed and integrated with respect to r|i so that 
finally there results 


U*. 

nb 


iL r £ 

211 Jo 


dT]2 


>/~T\ 2" 


*1] 


Gi(x,t) ) = - (b 2 — tj ) (q~ 


r 


dT)2 


(TI-TJ2 )</ Il2(l3 2 -Tl2)(a 2 - J n2) 


Using equation (29), we can now write 

r &2 dr)2 */(b 2 - j q) (Ti-a 2 ) 

G( x ^) = tr / . ■ . " ■ ■ + 

y ni('b 2 - j n2)(a 2 -q2) 


2p_ 



(31) 


Now using the equality 



G(x,t])dti 


it is possible to derive, after some manipulation, the relation 


. T = -|ib 


i k» 2 K - E + - (n - ^ b 2 k» 2 


( 32 ) 


where K and E are complete elliptical i ntegrals of the first and 
second kind, respectively, of modulus k =y 1 - k» 2 = 
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For the Kutta condition to apply at the trailing edge, equation (32) 
must reduce to the expression r = - y T>E> \i = - ay. in the interval 

— t<y<t (region II, fig. 5). This leads to the equation for the trail- 
ing edge in region II given by 


2kE 

rtk’ 2 


(33) 


Figure 7 shows a plot of the trailing edge derived from equation (33). 
The asymptotic slope of this edge is 45° and the asymptotic width of the 
wing (i.e., b-a as b— *•<») can be shown to be l/n. By means of this 
plot, it is possible to find the relation between aspect ratio and semi- 
span s, for wings of the type shown in figure 5. This was done numeri- 
cally and the results are shown for two different tip conditions in fig- 
ure 8. 


Begion Ili (x^s) .- Consider for a moment that the wing plan form 
ends with region II j that is, the wing has the form designated as type 
b in figure 8. The dovmwash in any plane behind the straight portion 
of the trailing edge of such a wing is determined by substituting equa- 
tions (31) and (28b) into equation (23) and integrating. Designating 
this value of w Q as w a , there results for t< |y| < s 


Wo, = yu ] 


and for 0 < |y| < t 


*00 = 3 ^ 


1 + 




where Ej. and F x are the incomplete 

modulus ki = t/s and argument y/t, and 


elliptic integrals with 
is the complete form of 



It should be noticed that behind the straight portion the flow 


is turning as if the wing were continued. Hence, any surface behind the 
straight part of the trailing edges of the type b wing will support 
zero loading, provided, of course, this surface is also rotating at the 
rate w = ny. Thus the solution for the type a wing is exactly the 
same as that for the type b wing with the added condition that the 
loading is zero over portions of the wing for which x>s. 


Besults for entire wing .- The results for the entire wing plan 
form can be summarized in the form of the total rolling moment and 
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subsequently the value of C7 . The rolling moment for both types of 

P 

wings is found from the expression 


L 


2pV c C 2 tan 2 0 



(35) 


Using the solution for T given by operating on equation (3l) , that is 
for 0 < | y | < t 


r = - ny 


and for t < y < s 


r 


r y2 . U_ / (s a -n)(n-t a ) J, 

2« J 2 J T)-qi J tii(s 2 -t)i) (ts-^i) 2 

s o 



*1 


in equation (35) and reversing the order of integration, the value of the 
moment can be obtained in coefficient form as 


3Cl p 
0 tan 0 


1 ^r 
4 0 tan 0 



5k! 2 - 1 
6s 


Ei 



(36) 


where again 0 has been. introduced as a mult iplier of spanwise and verti- 
cal lengths and where k x = t/s = Jl — kj/ 2 . Equation (36) is plotted 
in figure 9 for the plan-form-aspect— ratio relation given in figure 8. 

When k x is zero, equation (36) reduces to equation (21) for the triangu- 
lar wing. „ 


Arbitrary Plan Form 

If all the terms in equation (28a) are considered then equation (28b) 
becomes 


G*(m) 




I 

n = l 


n-g 

nBnTii 


(28c) 


and, using equation (29), equation (30) can be written 
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m 


H + 


!=?/ 

n=l o 


„2 2=5 

TJi 2 dT]i H(x) 1 


n^ii 


= -~ f 

2 it J 


G i( x >Tii)dti] 


(37) 


a® n-ni 


Now if 


m 


H(x) = 2n - gi \ nB n a 


2|i r . n— x p K n — i 


sn u du 


n=l 


then 


■^2 in 2 n— c 

Gi( x ,n) = |^(b 2 -n)( n -e 2 ) f V 

a 2 (^-^ 2 ) ^ 2 * t i-hi 


os) 


and, finally. 


r - ^ nB n a J (l+k 2 )sn n 1 xt-2k 2 sn n+1 u du — n V Bna* 1 — 

n L J <- — 1 


n=Q_ 


* a. 2 ^ m*: 


n=l 


m 


2n £ b k* + 


31 I BBn ’ >n 1 / 


K 


n — 1 
sn u du 


(39) 


n=l 


In order to satisfy the Kutta condition, r in equation (39) must equal 
m 

^ B n a n . Thus, the equation of the trailing edge becomes 

m 

^ nBna 11 ” 1 !!! (4o) 


n=l 


a = 


2k 

irk»2 


n=l 


and, using equation (35), the expression for rolling-moment coefficient 
becomes 


p % . v 


m 


n-i n — 1 , 2 
nB n k x s (VI 

■ 

n=l 


/ 2 t l-3ki‘ \ 

^ 1 I n +2 + 2 I xx) ~ 


3*k x » 

32 


(M) 
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where 



Values of the first six are given in Appendix B. 

Since the trailing edge must pass through the point a = 0 when 
x = 1, the value of Bi must always be 1. The’ other coefficients, 
however, axe at our disposal. These can be determined by specifying the 
points through which the trailing edge must pass and setting up a series 
of simultaneous equations using equation (4o). 

Two examples have been considered which represent a wing with very 
little taper (fig. 10(a)) and a wing which tapers to a point (fig. 10(b)). 
The trailing edge for the former of these wings, type c, was found by 
finding the values of B2 in equation (40) which made a = 0.4 when 
b = 1.2 and all the B* s with higher subscripts were zero. The trailing 
edge of the latter of the two wings, type d, was determined by finding B2, 
B 3 , B 4 , and B 5 such that the edge passed through four points which were 
on a straight line joining the root with the tip. Figure 11 shows the 
actual edge resulting for the type d wing together with the intended 
straight line and the points chosen. The slight scallop indicated in the 
trailing edge was found in the case of the triangular wing (i.e., when 
this method- was used to solve the delta— shaped wing and the results 
compared with the exact 'solution of the linearized equation) to give a 
3C Z 

value of ■ * ■ about 5 percent too low. However, this was considered 
p tan 0 

to be sufficiently accurate to establish the trend which is shown for 
this type of wing in figure 12. 


Ames Aeronautical Laboratory, 

Natio n a l Advisory Committee for Aeronautics, 

Moffett Field, Calif., July 22, 1949. 

APPENDIX A 

TABLE OF IMPORTANT S1MB0IS 

a, a(x) radius of cylindrical body in transformed space 

7 (in portion of report titled "Triangular Wing with 

Body" ) 

a, a(x), a*(y) equation of trailing edge of swept wing in transformed 

space (in portion of report titled "Swept-Back Wing," 
see equation (25)) 
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radius of cylindrical tody or equation of trailing edge 

reduced aspect ratio, */ |1-M 0 2 | times aspect ratio 
equation of wing leading edge in transformed space 
equation of wing leading edge 
root chord of wing 


rolling moment 
(2qS)(wing area) 

damping— in— roll coefficient 


c^c -1 

n, 

- B( P s/v 0 ) - 


damping— in— roll coefficient of wing alone 


Ei(x) 


G(x,tj) 


complete elliptic integral of the second kind with 
modulus k]_ 


incomplete elliptic integral 


1— k-, 2 t* 


function introduced in equation (24) 

complete elliptic integral of the first kind with modu- 
lus kx 


Fi(x) 


incomplete elliptic integral 


(l-k 1 2 t 2 )(l-t 2 ) 


rolling moment 


Mach number of free stream 
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P 

q 

r 

B 

s 

S 

sn(u) 

T, t 

v o 

U, V, V 


*t, n, wt 

x, y, z 
X, T, Z 

P 


q 


r 

n 

0 
n 

1 

p 

$ 


rate of roll, radians per second 
dynamic pressure ^ipVo 2 ^ 

radial length in polar coordinates 
ratio of "body diameter to wing span 
wing semispan in transformed space 
wing semispan 
.Jacobian elliptic function 

distance used in study of swept-back wing (see fig. 5) 
free— stream velocity 

transformed perturbation- velocity components in x,y,z 
directions, respectively 

radial component of velocity 

perturbation— velocity components in X,T,Z directions, 
respectively 

transformed Cartesian coordinates 
Cartesian coordinates in physical, plane 
</ | 1 - M 0 2 | ' 

loading coefficient, 

/ -pressure on lower surface - pressure on upper surface 

V q 

spanwise distribution of circulation 

y 2 

semiapex angle of wing 
pC 2 tan 2 0 

complex variable (y + iz) 
free-stream density 
perturbation potential 
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perturbation potential' in transformed coordinates 
jump in potential in the z = 0 plane 
polar angle 

APPENDIX B 

Exact expressions for lx, I 2 > . . . Ig. 


9 

A <P 0 

♦ 


lx = E 


_ 1 k» 2 . 1+k 

I 2 = - + t Zn 

2 4k 1-k 


I 3 = 3 ^- [ k» 2 K - (l-2k2) E] 


T = i _ i 
4 4 8 k2 


1 - 


l± 3 k 


Zn 

2 k 1 -k 


) 


Is = 


^5 2 ( 1 — 4k 2 ) ( 1 +k 2 ) 

5 + 15k2 


K-E l-4k 2 
k2 15k2 


I 6 = - 


( l+ 3 k 2 ) ( 3 — 5 k 2 )- ( 1 -k 2 ) ( l+ 2 k 2 + 5 k 4 ) 


48k< 


Zn 


32 k- 


1 +k 

1 -k 


Notice that vhen k=l, I n = — . 


Series expansions for lx, I 2 , . . . I g for k<l. 

Ii = — f 1 — — — ~k. 4 — k^- . . . ^ 

2 \ 4 64 2£6 5f / 


Is = 1 - - k 2 - — k 4 — k 6 - . . . 

3 15;’ 35 

I 3 = f ( 1 - | k2 - Ir- k 4 k 6 - 

3 4 \ 8 64' '■ 1024 


) 


Is = 


Ie = 


2 __ 

— k 2 - 

- — k 4 - 

— k 6 - . 

3 

15 

35 

315 

it 

( 3 -f 

k2- 

k4--i3 

16 

\ 4 

128 

512 

_8 

- k 2 

-i^k* 

_ i<L k s 

15 

35 

315 

693 


k 6 — . . 
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Figure I - Comparison of two theoretical results for a 
rolling triangular wing. 





Transformations 
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Figure 3 - Physical and transformed planes used in analysis of rolling wing with body. 
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Figure 4- Effect on damping in roll of a cylindrical 
body mounted on. a pointed wing with unswept 
trailing edge. 
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F igure 5 - Dimensions , transformations } and regions used in analysis of rolling 
swept- back wing. 
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Figure 7. - Graph of trailing edge y= a(x) for special swept- back wing. 
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Figure 9. -Variation of damping-in-roll parameter with reduced- aspect-ratio parameter 
for low values of reduced aspect ratio. 






